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ABSTRACT 

We investigate the effects of the magnetic dipole- dipole coupling and the 
electromagnetic radiation on the frequency evolution of gravitational waves 
from inspiralling binary neutron stars with magnetic dipole moments. This 
study is motivated by the discovery of the superstrongly magnetized neutron 
stars, i.e., magnetar. We derive the contributions of the magnetic fields to 
the accumulated cycles in gravitational waves as N mag ~ 6 x 10~ 3 (if/10 16 G) 2 , 
where H denotes the strength of the polar magnetic fields of each neutron star 
in the binary system. It is found that the effects of the magnetic fields will be 
negligible for the detection and the parameter estimation of gravitational waves, 
if the upper limit for magnetic fields of neutron stars are less than ~ 10 16 G, 
which is the maximum magnetic field observed in the soft gamma repeaters 
and the anomalous X-ray pulsars up to date. We also discuss the implications 
of electromagnetic radiation from the inspiralling binary neutron stars for the 
precursory X-ray emission prior to the gamma ray burst observed by the Ginga 
satellite. 

Subject headings: binaries: close — gamma rays: bursts — gravitation — stars: 
magnetic fields — stars: neutron — waves 

1. INTRODUCTION 

Direct detection of gravitational waves (GWs) is one of the most exciting challenges 
in the history of science. Long baseline interferometers for detection of GWs such as 
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TAMA300 (Kuroda et al. 1997), GEO600 (Hough 1992), VIRGO (Bradaschia et al. 1990), 
and LIGO (Abramovici et al. 1992) will be in operation within five years. One of the 
most promising sources of GWs for such detectors is the inspiralling binary neutron stars 
(BNSs), since we may expect several coalescing events per year within 200Mpc (Phinney 
1991, Narayan, Piran, & Shemi 1991, van den Heuvel & Lorimer 1996). 

As the orbital radius of BNSs decays due to gravitational radiation reaction, the 
frequency of GWs sweeps upward in detector's sensitive bandwidth from ~ 10Hz to 
~ 1000Hz. In the early inspiralling phase of BNSs, each neutron star (NS) can be treated 
as a point particle and the post-Newtonian (PN) expansion will converge (Cutler et al. 
1993) since the orbital separation of NSs is much larger than the NS's radius and the 
orbital velocity is much smaller than the velocity of light. This means that the theoretical 
templates of GWs in the inspiralling phase can be calculated with high accuracy by the 
PN approximation of general relativity using only several parameters: each mass and spin 
of NSs and the initial orbital elements. By cross correlating the observed noisy signals 
with the theoretical templates, the binary parameters, such as masses, spins (Cutler et 
al. 1993, Kidder, Will, & Wiseman 1993, Cutler & Flanagan 1994, Poisson & Will 1995), 
and cosmological distances (Schutz 1986), are determined. The number of rotation of 
BNSs is about 16000 in the detector's sensitive bandwidth. Therefore, the quite accurate 
theoretical templates are needed in order to extract physical information about BNSs 
from GWs since any effect that causes only one cycle ambiguity over 16000 accumulated 
cycles in the theoretical templates will reduce the signal-to- noise ratio (SNR). For the 
inspiralling BNSs in the sensitive bandwidth, with v 2 ~ m/r (hereafter G — c — 1) typically 
around 10~ 2 , the correction of 1/16000 ~ 10 -4 corresponds to second PN (2PN) order, 
(m/r) 2 ~ 10~ 4 (Blanchet et al. 1995). Many efforts are devoted to calculate higher-order 
PN corrections to theoretical templates (e.g., Blanchet 1996, Jaranowski & Schafer 1998a, 
1998b, Damour, Jaranowski & Schafer 1999, Tagoshi & Nakamura 1994, Tagoshi & Sasaki 
1994, Poisson 1995). However, these studies pay attention only to the gravitational effects 
on the theoretical templates, and it has not been studied how large corrections to the 
theoretical templates are caused by the electromagnetic effects, i.e., the magnetic fields of 
NSs, as far as we know.Q 

NSs observed as radio pulsars are believed to have strong magnetic fields, typically 
~ 10 12 G, assuming that the spin-down of pulsars is due to magnetic dipole radiation (e.g. 
Taylor, Manchester, & Lyne 1993). In addition, it begins to be recognized recently that 
NSs with superstrong magnetic fields ^ 10 14 G really exist (see below). We can crudely 
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estimate the correction to the waveform due to magnetic dipole fields of NSs by comparing 
the gravitational force F G ~ mim 2 /r 2 and the magnetic force F M ~ 3/ii/i 2 /V 4 between 
BNSs of masses, m\ and m 2 , with magnetic dipole moments, Hi and /i 2 , as 

E*L rr 4 f r ^ 2 ( El \ ( H * 

Fq "' X ^6(mi + m 2 )y V2 x 10 16 Gj V2 x 10 16 G 



x 



/ Ri \ 3 ( R2 \ 3 / 1.4M \ / 1.4M Q \ / 2.8M \ 2 

U0 6 cm/ U0 6 cm/ V mi / V m 2 )\m 1 +m 2 ) ' U 



where if p = 2/i p /Rp {p — 1,2) are the magnetic fields at the pole of the star and FL P are 
radii of NSs. This corresponds to the 2PN order correction. Therefore, the magnetic fields 
of order ~ 10 16 G might cause about one rotation error. Note that the r-dependence of the 
magnetic correction is the same as that of 2PN order, t> 4 oc r~ 2 , so that this argument is 
independent of the value of the separation r. 

Theoretically, in a new born NS, such superstrong magnetic fields ~ 10 16 G can 
be generated if the initial spin is in millisecond range since the conditions for helical 
dynamo action are met during the first few seconds after gravitational collapse (Duncan 
& Thompson 1992, Thompson & Duncan 1993). Observationally, such a superstrongly 
magnetized NS, or "magnetar", may be found as the soft gamma repeaters (SGRs) and the 
anomalous X-ray pulsars (AXPs) (Duncan & Thompson 1992, Thompson & Duncan 1996). 
The dipolar magnetic fields of SGRs are estimated as ~ 10 15 G by using the measured 
spin periods with spin-down rates for SGR1900+14 (Kouveliotou et al. 1999, Woods et 
al. 1999b) and SGR1806-20 (Kouveliotou et al. 1998), and with the peak luminosity for 
SGR1627-41 (Woods et al. 1999a). APXs that have measured the spin-down rate (Israel et 
al. 1999, Mereghetti, Israel, & Stella 1998, Baykal et al. 1998, Vasisht & Gotthelf 1997, 
Haberl et al. 1997) can be considered to have magnetic fields of 10 14 -10 15 G. Although 
there are some other models for SGRs and AXPs, these objects are best understood within 
the framework of magnetar (see, e.g., Thompson & Duncan 1995, Kouveliotou et al. 1998, 
Vasisht & Gotthelf 1997). Therefore, from both theoretical and observational results, it 
may be possible that a NS has superstrong magnetic fields ~ 10 16 G. 

In this paper we investigate the effects of the magnetic dipole fields of NSs on the 
frequency evolution of GWs from the inspiralling BNSs, since the magnetic fields of 
order ~ 10 16 G might cause about one rotation error. Throughout the paper, we use the 
relation H = 2fi/R 3 to connect the magnetic moment fi to the magnetic field at the 
magnetic pole H, and R = 10 6 cm as the radius of a NS. For later convenience, note that 
/I = 1.4 x 10 9 (F/10 16 G) cm 2 in units of G = c = 1. 
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2. EQUATIONS OF MOTION 

We consider a binary system of two compact bodies of masses, m\ and m 2 , with 
magnetic dipole moments, ^ and ^i 2 , respectively. Since we pay particular attention to 
the effects of magnetic fields, we treat the orbital motion of BNSs in Newtonian gravity. 
Although a spherical symmetry of the stellar configuration is, in general, incompatible with 
the presence of magnetic fields (Chandrasekhar 1981), we ignore quadrupole effects which 
are caused by magnetic fields for a moment (see §[| and §|B|). We also neglect tidal effects 
which are expected to be small until pre-merging phase of BNSs (Bildsten and Cutler 1992). 

By eliminating the motion of the center of mass of BNSs and setting the origin of the 
coordinate frame at the center of mass of BNSs, the effective one-body equations of motion 
can be derived from a Lagrangian, 

£ = l V mv 2 + + c DD , (2) 
2 r 

where 

Cdd = Mi • H 2 = 4 {3(™ ■ ■ M2) - Mi ■ M 2 } • ( 3 ) 

Here m = mi + m 2 , rj = mim 2 /m 2 , r = \x\, x = X\ — x 2 , fi = x/r, v = x, and 

H 2 = {3(n • /x 2 )n — n 2 } / r3 i s the magnetic field at x± produced by the magnetic moment 

H 2 - By using the Euler-Lagrange equations, we obtain the equations of motion as 

in 

a = -—h + a DD , (4) 

where a = x and 

3 

a D D = 7 {Oi • M 2 )™ + (n ■ M 2 )Mi + (n ■ Mi)m 2 - 5 (™ ■ Mi)(™ ■ M 2 )™} • ( 5 ) 

From equation (|2]), the energy of this system is given by 

1 2 7/m 2 

= -rjmv h E DD , (6) 

2 r 

where 

£dd = — $ {Hn • • th) - Pi • Ih) • ( 7 ) 

The total angular momentum can be defined as L = + S where = rpn(xXv) 
is the Newtonian orbital angular momentum and S = Si + S 2 is the total spin angular 
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momentum. We can show explicitly E = L = with the equations of motion (^) and the 
evolution equations of the spins] 2 ] 

S p = (J, p xH q = — {3(ii- fj, q )(fj, p xh) - n p XfJ, q }. (p,q = 1,2) (8) 

Assuming NSs as spherical compact bodies, the spin angular velocities O p are related to 
the spins S p as S p = I p £l p where I p is the principle moment of inertia of the bodies. Since 
the magnetic moments evolve as fi p = Q P X fi p = S p Xfi p /I p , the angular velocities of the 
magnetic moments £l p will be of order £l p ~ (fXi^/mR 2 ^) 1 ^ 2 from dimensional analysis 
with equation (K). Note also that the orbital angular velocity w is of order w ~ (m/r) 3,/2 /m, 
and the orbital inspiral rate w ins = (dE/dt)cw/ E is of order w ins ~ (m/r) 4 /m. 



3. GRAVITATIONAL WAVES AND ELECTROMAGNETIC WAVES 



As a first step, we use the quadrupole formula to derive the rate of energy loss from 
a binary system due to GWs (e.g. Thorne 1980). The symmetric, trace- free parts of the 
quadrupole moments of this system are given by fy = r\m{xiXj — ^r 2 5ij). Taking third 
time derivatives of these quadrupole moments, we obtain the energy loss rate from the 
quadrupole formula as 



'dE\ 
,dt J Gw 



"5 V y ^' y 
8 ?7 2 m 4 I 



15 r 4 
1 



I2v 2 - llf 2 



rjm 2 r 2 



6(-12v 2 + 13^)0! • // 2 ) + 12(21t» 2 - 34r 2 )(n • /x x )(n ■ /x 2 ) 



-36(« • fJ,x)(v ■ /j, 2 ) + 87f{(n • ^)(v ■ /x 2 ) + (v ■ ^{n ■ /x 2 )} 



(9) 



where we have used the equations of motion (f|) and assumed fMi^/m r <^ 1. 

On the other hand, electromagnetic (EM) waves are also emitted from this binary 
system since the magnetic moments are moving. Using the linearity in the EM fields, the 
radiation fields B™ d in equation (|A4j) for this binary system are given by 



B 



rad 



D 



(d- p eff ){{d-a)d-a} 



(10) 



2 Note that there are no spin-orbit and spin-spin interactions since we consider Newtonian gravity. The 
spins are generated by the torque due to the magnetic dipole-dipole interaction. 
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where 



Veff = — -m x fx 2 ). 
11 m 



Therefore, since the radiated power is given by equation 
to the EM radiation is calculated as 



(11) 

the rate of energy loss due 



'dE\ 



EM 



2 m 2 
15 r 6 



2/Ue// {v 2 - 6r(n • v) + 9f 2 } - {fi eff • (v - 3rn)}' 



(12) 



where we have substituted fx e ^ into fx in equation (|A6|) . 

Note that the assumption of the constant magnetic moments in equations (|D|), flTU| ) and 
( |A4l ) is valid since the angular velocities of the magnetic moments Vt p ~ (fii^/mFPr 3 ) 1 / 2 are 
much smaller than the orbital angular velocities w ~ (m/r) 3//2 /m when (HiH 2 ) l l 2 <^ 10 18 G. 



4. INSPIRAL OF CIRCULAR ORBITS 



For calculational simplicity we assume that the orbital motion of BNSs has decayed 
to be circular apart from the adiabatic inspiral (Peters Sz Mathews 1963, Peters 1964). In 
general, circular orbit solutions of equation (|4]) do not exit unless the magnetic moments 
are aligned perpendicular to the orbital plane. However, since the angular velocities of 
the magnetic moments fi p ~ (ni^/rnl&r^) 1 / 2 are much smaller than the orbital angular 
velocity w ~ (m/r) 3 / 2 /m when (HiH 2 ) 1 ^ 2 10 18 G, we can regard the magnetic moment 
vectors and L as time-independent ones over an orbit where L is a unit vector orthogonal 
to the orbital plane. Then, after taking an average of the magnetic term in the acceleration 
([]), we can obtain orbits of constant separation, r = r = n-v — 0, w — v/r, Ln = r]mr 2 wL 
(see similar discussions on spin precessions, Kidder, Will, & Wiseman 1993, Kidder 1995). 



From the equations of motion for circular orbits, n • a 
orbital angular velocity as 



r — rw 2 , we can calculate the 



m 



w 



1 + 



2r]m 4 



{Mi • - 3(Xr • fJt^L ■ fx 2 )} 



where we have used the orbit-averaged relation, 

(n ■ fx x ){n ■ fx 2 ) = - {[x± ■ fx 2 -{L- fx^L ■ fx 2 )] . 



(13) 



(14) 



The total energy and the energy loss rate for a circular orbit, averaged over an orbit, 
can be obtained as 



E 



n 



m 
2r 



1 

2i]m 4 



m 



{mi ' M2 - 3(X ■ fx x ){L ■ fx 2 )] 



(15) 
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dE 



'dE s 
dt 



+ 



GW 



'dE s 
~dt . 



32 2 / m \ ; 



5 ^ V r 



7 



9 / m x 2 



y / fft \ r 1 

1 + 2^l7i l^i-^-s^-^)^-^)} 



+ 



/// 



(16) 



— r/m 5/3 w 11/3 {l + 



)4/3} 



(17) 



96r/ 2 m 4 \ r 
by using equations @, @, (|P2|) and (|13|). 

Combining equations (|T3"D, flI5| ) and QTfiP, we can express the change rate of the orbital 
angular velocity zi; as a function of w, 

96 

y 

where 

"mag = ^ {Ml • M 2 - 3(X • ' M 2 )} + + (Me// • X) 2 } • (18) 

By using equation (|T7|), we calculate the accumulated number of GW cycles N = J(f/f)df, 
where / = w/n is the frequency of the quadrupolar waves. For calculational simplicity we 
assume that the two magnetic moments are aligned parallel to the orbital axis of the binary 
system. In this case, a mag becomes a constant value. Then, the accumulated number is 
integrated as 

(19) 



N = N + N 

ly ly grav ~ 1 ' mag ; 



where N grav and N mag denote the contributions from the Newtonian gravity term and the 
magnetic term respectively, and are expressed as 

-i fmax 



N. 



grav 



N, 



32nr] 
5 



[irmf) 



-5/3 



(20) 



fmin 
.fn 



mag 



-a 



32tt?7 



mag 



[irmf) 



-1/3 



(21) 



Here f max is the exit frequency and f min is the entering one of the detector's bandwidth. 



5. DISCUSSION 

Using 10Hz as the entering frequency and 1000Hz as the exit one, we obtain the 
contribution to the total number of GW cycles from the magnetic term as 

, / //. \ / H 2 \ / m x " J:!/:{ 



N <~ = - 5 - 9 x 10 Iio^gJ Kwk) [jtmrj ' ^ 
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where we assume n l ■ fi 2 < 0, /x 1 || L, \x 2 || £,0 and m\ = m 2 . Note that the contribution 
of the EM radiation reaction (the second term in equation flI8"D) is much less than that of 
the dipole-dipole interaction (the first term in equation flUf)). The maximum magnetic 
field allowed by the scalar virial theorem is ~ 10 18 G for NSs (Chandrasekhar 1981, see 
also Bocquet et al. 1995). If NSs in the inspiralling BNSs have such magnetic fields 
10 17 G <; H <^ 10 18 G, the effects of the magnetic fields can change more than one cycle in 
the accumulated cycles. However, if we consider that the maximum value of the observed 
fields ~ 10 16 G is the upper limit for the magnetic fields of NSs, the magnetic term will 
make negligible contributions to the accumulated phase, contrary to the crude estimate in 
§|Tj and equation ([I]). Consequently, the magnetic fields of NSs will not present difficulties 
for the detection of GWs from BNSs, if the upper limit for the magnetic fields of NSs is less 
than ~ 10 16 G. 



The magnetic term in equation (|T7|) has the same dependence on the angular velocity 
w as 2PN terms. We can see this dependence from the 2PN expression for the frequency 
sweep (Blanchet et al. 1995), 



w 



r]m 5/3 w 11/3 jl - + ^-r^J (mw) 2/3 + (47T - (3)(mw) 



9_6_5/3...n/3 fi ^43 , 11 
5 

/34103 13661 59 2 \ . , 4/ „l 

+ 1 77 + — v + o (mw 4/3 \ , 23 

V 18144 2016 18 / J V ; 

where the spin-orbit parameter is (3 = ^S p (113mp/m 2 + 75r])L ■ x p , the spin-spin parameter 
is er = (77/48) (—247% 1 • x 2 + 721L ■ X\L ■ X2) an d X P — S p /m 2 . As we can see from equations 
( P7|) and (PB|), we cannot distinguish the magnetic term a mag from the spin-spin term cr at 
2PN order. Only the sum of the magnetic term and the spin-spin term can be deduced 
from the frequency evolution of 2PN order. This degeneracy might be broken by examining 
the wave-form modulations caused by the spin-induced precession of the orbit (Apostolatos 
et al. 1994, Kidder 1995). However, the analysis of Poisson & Will (1995) shows that 
the measurement error on the spin-spin term is Act ~ 17.3 assuming that the SNR is 10, 
(3 = 0, a = and there is no prior information.^ On the other hand, the contribution of the 



3 This set of configurations, fi ± ■ fi 2 < 0, fi 1 L and /x 2 || L, is the most stable one, 
i.e., the EM interaction energy in equation (0) becomes the minimum value. However, note that 
the alignment rate of the magnetic moments w a u is smaller than the orbital inspiral rate Wi ns ~ 
(m/r) 4 /m when H < 10 16 (r/10 n cm) ' G, since the alignment rate w a n can be estimated as w a n ~ 
(energy loss rate due to magnetic dipole radiation)/ Eod ~ (M^p)/^ 2 / 7 " 3 ) ~ /x 4 /r 3 m 2 i? 4 when \i\ ~ ^2 ~ 

4 Because of some simplifying assumptions in Poisson & Will (1995), the true measurement error is still 
uncertain. However, this will not affect the discussion since the true measurement error will be within a factor 
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magnetic term is estimated as 



a mag = 2.9 x 1(T 3 f-gU f-gf-) ( - - ™ - ) , (24) 



Vio 16 g; uo 16 g/ v 2 x 1-4M 0/ 

where we assume /x 1 • /x 2 < 0, /x x || L, n 2 II an d m i = m 2- We can see from this equation 
that the contribution of the magnetic fields of NSs is much smaller than the measurement 
error on the spin-spin term. Therefore, the effects of the magnetic fields of NSs are also 
negligible for parameter estimation with moderate SNR if we consider that the maximum 
value of the observed fields ~ 10 16 G is the upper limit for the magnetic fields of NSs. 

Conservatively, there are some other reasons to consider that the magnetic fields of 
inspiralling NSs are smaller than ~ 10 16 G. There are three observed BNSs in our Galaxy 
and nearby globular cluster that will merge in less than a Hubble time: PSR B1913+16 
(Taylor & Weisberg 1989), PSR B1534+12 (Wolszczan 1991, Stairs et al. 1998) and PSR 
B2127+11C (Prince et al. 1991). From the spin-down rate, the magnetic field strength is 
estimated as of order 10 10 G for all pulsars in these binary systems. If such a value is typical 
of the magnetic field strength^ the magnetic terms will be negligible. Moreover, if the decay 
time of the magnetic fields is shorter than the coalescence time (Heyl and Kulkarni 1998, 
Thompson and Duncan 1996, Goldreich & Reisenegger 1992, Shalybkov and Urpin 1995), 
of course, the magnetic fields will not be concerned, although the magnetic field evolution 
of isolated NSs is an unresolved issue. It may be also difficult for BNSs including magnetar 
to be formed because of a large recoil velocity (Duncan & Thompson 1992). 

In this paper, we regard NSs as spherical compact bodies. However, when we consider 
NSs as extended bodies, we have to take into account of the quadrupole effects induced by 
magnetic fields. The magnetic fields are a source of non-hydrostatic stress in the interiors 
of NSs. A magnetic dipole moment /x would give rise to moment differences of order 



I C — la R H 



(25) 



where I c and I a refer to the moments of inertia about the dipole axis and about an axis in 
the magnetic equator. Then, the gravitational potential 7?m 2 /r between NSs is modified 
by an amount of order ml a e/r 3 ~ /i 2 /r 3 . This is the same order as the EM interaction 
term in equation (E3) when \ix ~ Therefore, an accurate ellipticity e in equation (EH) is 



of the order of 2 (Balasubramanian, Sathyaprakash, & Dhurandhar 1996, Balasubramanian & Dhurandhar 
1998, Nicholson & Vecchio 1998). 

However, note that no radio pulsars have magnetic fields above ~ 10 G, and hence there may be a 
selection bias (Baring & Harding 1998). 
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needed to determine the magnetic effects within a factor. (Bocquet et al. 1995, Bonazzola 
& Gourgoulhon 1996, Konno, Obata, & Kojima 1999). Note that there are such quadrupole 
effects even if only one companion has magnetic moment. 

Although the effects of the magnetic fields of NSs will be negligible for observations 
of GWs, they might be concerned with gamma ray bursts (GRBs). The BNS 
merger is one of models of GRB sources, (see, e.g., Piran 1999 for a review). If 
a NS in the binary system has strong magnetic fields, the total energy emitted by 
EM waves until coalescence can be estimated from equations flT3|), (|i~6|) and (|17|) as 
-J(dE/dt) EM dt ~ Tr 2 ntff(fLx ~ fLJ/^Vm ~ 10 46 (tf/10 16 G) 2 ( Wl0 3 Hz) 2 ergs. 
This energy will be radiated at very low frequency ~ 10 3 Hz, which is difficult to be 
observed by the present radio telescope. Furthermore such radiation cannot propagate 
a plasma if an electron density is larger than ~ 0.01cm -3 since the plasma frequency 
is larger than the radiation frequency (e.g. Spitzer 1962). However, this energy may 
be converted to the thermal energy of the surrounding plasma efficiently if the electron 
density n e is sufficiently high since the electron-electron relaxation time is about 
t re i ~ l(ra e /10 11 cm~ 3 )~ 1 (/cT e /2keV) 3//2 s where T e is the electron temperature (e.g. Spitzer 
1962). Then, this thermal radiation might explain the precursory X-ray emission ~ 10 s 
before the onset of the GRB observed by the Ginga satellite, in which the total energy of 
the X-ray precursor emission is estimated to be about ~ 10 46 (<i/100Mpc) 2 ergs (Murakami 
et al. 1991). Even though the strong magnetic fields are not relevant to the X-ray precursor 
in GRB, the EM radiation can be the EM signature of the coalescing BNSs. Therefore 
it is an interesting future problem to investigate the conversion of the low frequency EM 
radiation to the higher frequency one. 
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A. ELECTROMAGNETIC RADIATION FROM A MOVING MAGNETIC 

MOMENT 

In this section we review the radiation from a moving magnetic dipole moment. We 
consider a particle with only a magnetic dipole moment fi' in its rest frame K' . A moving 
magnetic dipole moment with velocity v = x relative to an observer frame K also has an 
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associated electric dipole moment. The apparent electric dipole moment is 

p = vxp, (Al) 

where p = p! — -^{v • p')v is the magnetic moment observed in K and 7 = (1 — t> 2 ) -1 / 2 
(Jackson 1998). Therefore, in the observer frame K the magnetization density M and 
electric polarization density P are given by 

M(t,z) = p(t)6[z - x(t)), 

P(t,z) = p(t)5[z-x(t)\, (A2) 

where p is given by equation (|A1|). Recalling that the moving magnetic moment is 
equivalent to a current J = VxM + P (Jackson 1998), we can calculate the electric and 
magnetic fields. For our purpose, it is sufficient to obtain the radiative parts which fall 
off as the inverse of the distance D~ l . The radiation field B rad of this moving magnetic 
moment is given by (e.g. Heras 1994) 

d 3dx(dxp — p)(d- a) 2 3dx (dx/i - p){d ■ a) 



B ™ (t ' z) ~ DW 5 ' DW 4 

dx(dxp, - p)(d ■ d) ( dx(dxfl-p) o 

(A3) 



DW 4 DW 3 

where D(t) = z - x(t), d = D(t)/\D(t)\ = D(t)/D, W = 1 - v • d, a = v and the 
right-hand side of this equation is evaluated at the retarded time t', i.e., t' + D(t') = t. 
When we assume that the magnetic dipole moment vector p,(t) is a constant one, the above 
equation (|A3|) can be calculated as 



rtrad 

B ° : D 



i |dx(dx^)(d ■ a) - dxpj 
= i(d-/x){(d-d)d-d}, (A4) 

up to the leading order term of v and \x\/D. Note that the radiation field B™ d becomes 
the same equation as (|A4|) up to leading order even if we use p! instead of p. The power 



radiated per unit solid angle is given by (e.g. Landau & Lifshitz 1975) 

The total instantaneous power is obtained by integrating equation ( |A5| ) over all solid angle 

as 

P = -{2 M 2 |d| 2 -(/x-d) 2 }. (A6) 
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B. THE CONTRIBUTION OF ELECTROMAGNETIC FIELDS TO THE 

MASS 

Thus far we have ignored the contribution of the EM fields to the mass of a compact 
body in a binary system. In this section, we estimate the correction to the mass by the 
EM fields. We have implicitly defined the mass m as that of the isolated spherical body, 
i.e., the orbital separation of the binary system is infinity. Therefore, the mass m includes 
the self-energy of the EM fields. The self-energy m ma g of the magnetic fields outside the 
compact body with a magnetic dipole moment /j, is obtained by 



f — ^ > - — f 3/i 2 COS 2 & + /i 2 ^3 , _ 2 

m mag -J ^dx -^y dx -fi 



3r /3 



,2 



R 



(Bl) 



3i? 3 



where R is the radius of the compact body. (The effects of the magnetic fields inside the 
compact body are discussed in §|[) 

Next, we consider the gravitational interaction between the above compact bodies 
with magnetic moments. Here we note that at a finite separation the gravitational field 
of a mass with magnetic dipole moment [i is different from that of a point mass without 
magnetic dipole moment even if the masses are the same value since EM fields have an 
extent. Therefore, the mass m in the gravitational potential term rjm 2 /r in equation (|2|) 
suffers a small correction. This correction can be estimated by evaluating the "gravitational 
potential" produced by the energy of the EM fields of the magnetic dipole moment \i. The 
"gravitational potential" at x is obtained by 



*) = - o ,: _„ <*v = r „ <*v- (B2) 



H 2 „ . ft 2 r 3 cos 2 9' + 1 ., 

-d A x = -— / — —. — 

87r|a; — x \ 8% J r' b \x 

Here we expand l/\x — x'\ by the spherical harmonics Yi m (8, ip), 

1 oo i i r i 

4^EE iTT-feW^v), (B3) 



tJC in i 2 / I 1 

i i (=0 m=—l 



where r < (r > ) is the smaller (larger) one of \x\ and \x'\. Noting the orthonormality 
of the spherical harmonics, f Y l f m ,(9,{p)Y[ m ( y 9,{p)dQ = 5i>i5 m r m , and a relation, 
3cos 2 fl' + 1 = 2yf^/5Y 2 o(9',(p') + 2V^Y 00 {6', if'), the equation QB3|) is calculated 

as 



<P(x) = -- 
r 



jj, 2 /j 2 (3cos6 — 1) /j 2 cos 2 6 



3R 3 lORr 2 4r 3 



(B4) 



The first term in the bracket on the right-hand side of the above equation comes from 
the total self-energy of the EM fields in equation ( |B1| ), and the last two terms are the 
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corrections due to the extent of the EM fields. The order of the correction to the mass 
is estimated as 5m/m ~ (/i 2 /Rr 2 )/m ~ fi 2 /mRr 2 . On the other hand, the correction 
due to the dipole-dipole interaction is of order /x^/mV 2 from equation Therefore, 
the correction due to the extent of the magnetic fields are smaller than that due to the 
dipole-dipole interaction when /ii ~ fi2- We have also confirmed that the contribution of 
the EM fields to the quadrupole moments is of order n 2 /mRr 2 . 
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